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Abstract 

We present a strong field theory of matter wave splitting in the pres- 
ence of various gravitational, inertial and trapping potentials. The effect 
of these potentials on the resonance condition (between the splitting po- 
tential and the considered effective two-level system) and on the atomic 
Borrmann effect is investigated in detail. The dispersive structuring of 
an incident atomic wave packet - due to such generalized beam splitters - 
is studied and modeled, and several important dynamical features of the 
solutions are detailed (generalized Rabi oscillations, velocity selection, 
anomalous dispersion, generalized Borrmann effect and anomalous gravi- 
tational bending). Finally, we show how to express this triple interaction 
"matter - splitting potential - gravito-inertial and trapping potentials" 
as an equivalent instantaneous interaction which turns out to be a very 
efficient tool for the modeling of atom interferometers. 

PACS number(s): 03.75.-b, 32.80.-t, 33.80.-b, 39.25.+k, 42.50.-p 

1 Introduction 



The matter wave beam splitters are nowadays the cornerstone of a wide range 
of experiments, from atomic clocks and gravito-inertial sensors to laser cooling 
and ultracold atoms characterization, quantum computing and cavity QED ex- 
periments, atom lithography and chemical reaction dynamics, detection of tiny 
effects of General Relativity and test of fundamental theories, measurement of 
atom surface interactions... 



In view of the recent progress in non-dissipative atom optics (coherent beam 
splitters, mirrors, lenses...) as well as in dissipative atom optics (slowing, trap- 
ping and cooling of atoms and molecules), it is needed to deepen our compre- 
hension of light-matter interactions in the presence of other external potentials, 
like gravito-inertial or trapping potentials. In particular, the precision and sta- 
bility of atom interferometers are now so outstanding 1 that it is necessary 
to go beyond the former modeling of their main component, namely the beam 
splitters. 

In fact, the concept of atomic beam splitter is not confined to light-matter 
interactions, and can be extended to any interaction process between matter 
waves. It is thus possible to write the action of such atom optical elements 
as a S matrix between the incident and diffused matter waves, where the S 
matrix depends mainly on the "splitting" potential, which can be material (slits, 
periodic microstructures...) or electromagnetic (magnetic or electric static fields, 
laser fields...). This S matrix description is particularly useful for the modeling 
of atom interferometers, and more generally for any set up having a succession 
of such beam splitters E3 0] . 

However, for a long time, the precision of atom optics experiments has re- 
mained low enough not to require an accurate study of matter wave beam split- 
ters. Thus, in the most common simplified modeling of these elements, only the 
following effects were considered: 1) the splitting of an incident atomic wave 
packet into several wave packets, 2) among them one was equal to the incident 
wave packet, up to a change of amplitude, 3) and where the others could differ 
from the incident wave packet in their central momentum, internal state, ampli- 
tude and phase. However, this practical modeling - sometimes called "infinitely 
thin" because it amounts to neglecting the duration of the interaction - does 
not take into account several important effects, like the dispersive structuring 
of the incident wave packet (velocity selection and sidebands, Borrmann effect, 
anomalous dispersion...), the time and space dependency of the splitting poten- 
tial, the effect of relaxation processes, or the effect of other external fields during 
the splitting (time-dependent gravito-inertial effects, trapping potentials...). 

During the last past two decades, several authors studied some of these 
problems, namely: the effect of a non-trivial time dependency of the splitting- 
potential (for running and standing laser beam splitters) [5JEJ[3|S1|2|; the atomic 
Borrmann effect and anomalous dispersion effect without any other external 
potential ^3 or with a constant and uniform acceleration (WKB solution 
|12|h the atomic splitting in a constant and uniform acceleration (exact solution 
in the temporal case and WKB solution in the spatial case) a common 
modeling for both spatial and temporal beam splitters to the first order in the 
splitting potential (weak field theory) |14|... 

In the light of what happened in neutron optics, where the beam splitters 
modeling proved to be crucial to understand properly the origin of the inter- 
ferometer phase shifts ^3 E| j h appears to be necessary to go beyond these 
studies, so as to provide a comprehensive modeling of the true action of a matter 
wave beam splitter (strong fields theory for all the involved external fields). 

This paper is organized as follows. First, we give some details on our frame- 
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work and explain how to put in equation the problem of the triple interaction 
"matter - splitting potential - other external fields" . Then, we detail how to 
transform the obtained equation in a simpler one thanks to unitary transforma- 
tions (interaction picture) and passage into the rotating frames. We expound 
then how to solve this equation (analytically or numerically, with different de- 
velopments or relevant approximations) , and we go back to the initial represen- 
tation to explain how to write the effect of such beam splitters as an effective 
instantaneous interaction (generalized ttt scheme). Finally, we study the atomic 
Borrmann effect and other anomalous dispersive properties and model them in 
the general framework detailed in the second part. 



2 Framework and approximations 

2.1 General framework 

The matter wave beam splitters we consider in this paper consist of multi- 
level atomic systems subject to an interaction potential which couples the levels 
together. This interaction is usually made in the presence of other external 
potentials and miscellaneous relaxation processes. 

In fact, these "atomic systems" can refer to atoms (neutral or not) as well as 
molecules, and more generally to any quantity of matter which can be coherently 
manipulated. 

Furthermore, the atomic levels are not restricted to internal atomic levels, 
but more generally refer to energy- momentum states (i.e. eigenstates of both the 
internal and kinetic Hamiltonians) . The transitions can occur between internal 
states only (spectroscopy without Doppler effect for example), external states 
only (diffraction in Kapitza-Dirac and Bragg regimes, optical Stern-Gerlach 
effect, magnetic atom mirror...), or between entangled states, where the entan- 
glement may be between the internal and external states (stimulated Raman 
transitions for example) or between the previous energy-momentum states and 
the eigenstates of the interaction potential (Fock states of the quantized elec- 
tromagnetic field for example). 

There are many techniques to coherently split a matter wave, and each of 
them corresponds to a particular kind of splitting potential. Like the other atom 
optics elements, matter wave beam splitters use essentially two properties of 
atoms (or molecules): their wave property, and their interaction with external 
fields, electromagnetic or material. To date, the demonstrated matter wave 
beam splitters are based on: 

1. atom-matter interaction: front wave division (material slits) |17j . ampli- 
tude splitting (transmission material gratings) [TS|. reflection at crystalline 
surface JHj > quantum reflection [20] • • • 
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2. interaction with static magnetic or electric fields (for atoms having a dipo- 
lar electric or magnetic moment): longitudinal Stark effect |5T], transverse 
and longitudinal Stern-Gerlach effects magnetic mirrors [231; mirror 
for polar molecules [2], Y shaped magnetic guides |2SJ. . . 

3. resonant or quasi-resonant interaction with laser fields: reflection and 
diffraction by standing |26|. running |27U28j or evanescent [2H| laser waves 
(with spatial and/or temporal working), optical Stern-Gerlach effect |30) . 
stimulated Raman effect j2U and its derivatives (adiabatic transfer 31 , 
STIRAP [32], CHIRAP |33]...), magneto-optical beam splitters 0, X 
shaped dipolar guides 33 ■ ■ ■ 

In this paper, we will focus on this third kind of interaction, and more 
generally on the beam splitters for which the two-beam approximation is valid. 

As for the relaxation processes, they refer to all the processes which lead to a 
loss of coherence and/or a loss of atoms (spontaneous emission, inter- atomic col- 
lisions, absorption and interaction with the material microstructures. . . ). When 
they can not be neglected, the use of a density operator formalism is needed. 

Other external potentials may be present during the matter wave splitting: 
inertial and gravitational fields, trapping potential, van dcr Waals and Casimir 
potentials. . . In this paper, we consider all the time-dependent potentials which 
are at most quadratic in position and momentum. The corresponding Hamilto- 
nian is therefore: 



\ ^ JJ1 ^ —X, ^ > 

Hext = — p .(3 (t).p - — r .7 (t). r - r .a (t).p -mg (t).r + f (t) . p 

This includes the effect of non-uniform accelerations (~g (t) and 7 (i)), ro- 
tations (with an angular velocity fl (t) such that a (t) .~u := — fl (t) X it for 
any vector ~u), trapping potentials (— 7 (t)), non-zero curvature tensor (7 (t)), 

gravitational waves in Fermi's gauge (7 (t)) or Einstein's gauge (/3 (£)), and all 
the electromagnetic potentials which can be written as a development at most 
quadratic in position and momentum ("g*(i), 7 (t), f (tj). Furthermore, to 

keep an overall approach, the coefficients of H ext are time-dependent, and a, (3 
and 7 are expressed with non-diagonal 3x3 matrices. 



2.2 Approximations considered in this paper 

It is often possible to simplify this general framework and obtain an evolution 
equation between only two effective states by using some justified approxima- 
tions. 

First, the two-beam approximation is indeed valid when only two energy- 
momentum eigenstates are coupled. The coupling can be direct (for true two- 
level systems) or indirect (Raman transitions, spatial beam splitters in Bragg 
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regime). In fact, one can show that any N-photon transition (thanks to sev- 
eral running or standing laser waves) of a multilevel atom is equivalent to an 
effective 1-photon transition between two atomic levels if the other levels can 
be adiabatically eliminated |2)J f37\ . The effective photon may not be real. For 
example, in the Bragg regime, the wave vector of this effective photon is equal 
to Nh k and its frequency is equal to 0. Eventually, the spatial and temporal 
structure of the true laser beams appears only in the amplitude of the effective 
running laser beam. 

Second, the laser fields are considered as classical (coherent states of the 
quantized electromagnetic fields), but the calculations which follow are also 
valid for a transition between two dressed states |38j . 

Third, we suppose that the two atomic levels have a long lifetime and we 
neglect all the relaxation processes listed before. However, the instability of 
these levels, due to spontaneous emission, can be taken into account in an ap- 
proximative manner adding a non-Hermitian part to the atom-laser Hamitonian 



In addition, V em is chosen equal to the usual dipolar electric Hamiltonian 
(without spin) and we suppose that the other external fields are sufficiently 
weak to neglect their effect on the atomic levels and laser fields. 

Finally, the triple interaction "laser - matter - other external fields" can be 
written as a Schrodinger equation concerning two atomic states coupled by an 
effective running laser wave: 



where r op and r op are the position and momentum operators, and Hq the 
internal Hamiltonian (Eb > E a ): 



where u>q = [E^ — E a ) /h is the atomic transition frequency and 03 the usual 
third Pauli matrix. 

It is also possible to account for some relativistic effects by introducing two 
different masses [QQ]. For simplicity however, we will not take into account 
these relativitic corrections and use only one atomic mass in what follows. 



3 Interaction picture and rotating frames 

It is generally impossible to solve directly the equation Q (non-trivial time 
dependence of the right hand side, presence of two non-commuting operators 
r op and r op ), but it is possible to simplify it with the help of well chosen unitary 
transformations 0]. 



V em M 



ih—\V(t)) = (H + H ext (r^,p^,,t) + V em (r^,t))\^{t)) 



(1) 
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The main idea of this series of transformations is to eliminate progressively 
the different sources of evolution (internal and external) of the right hand side 
of Q). As each unitary transformation corresponds to a change of frame, we can 
see this succession of transformations as a succession of frame changes which 
aims at reaching the proper frame of the atom, or, at least, at reaching a "least 
movement frame" for the atom (from external as well as internal points of view) . 
In this especially suitable frame, it is easier to solve the evolution equation and 
several important pieces of information about the solution can be directly seen. 

First, let us go to the interaction picture with respect to Hq and H ext : 

\*(t)) = U 1 (t,t 1 )\<p 1 (t)) (2) 

with : 

U l {tM) = e-^ Ho - ( - t - tl) U ext (i,ii) 

and: 

Uext (Ml) = T ^exp ^-1 jf H ext {r^p,p^ P ,t') dt ^j 

where T is the chronological Dyson operator and t\ an arbitrary time (dif- 
ferent from t by definition). 
The equation becomes: 



ihj t \yi (t)) = e i*M*-*i V em (1C P (t, h) , t) e -**°.(t-ti) \^ (*)> (3) 

where R op is defined as: 

Rop (t, tx) = U ext (t, ti) r op U ext (t, ti) 

The external Hamiltonian H ext being at most quadratic in position and 
momentum, R op depends linearly on v p and Pop'' 

Rop (t, h) = A (t, t\) r^, + B (t, t\) Po^,/m + £ (t, t\) 

and is simply obtained through the classical solution of the Hamilton's equa- 
tions. One can show that the matrices A and B depend on the quadratic terms 

of H ext (a (t), (3 (t) and 7 (i)) only, whereas ^ (Mi) a l so depends on its linear 
terms (if (t) and / (t)). These matrices are in fact the well known ABCD ma- 
trices, usually used in Gaussian optics, and introduced recently in atom optics 

[2313]. 

As far as V em is concerned, it may have diagonal terms (AC stark shifts, 
slowly varying in space and time) . However, in a first approach, we can take the 
latter constant and eliminate their common part by a unitary transformation. 
Finally, V em can be taken as purely anti-diagonal: 
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with: 



1 

1 



V (r, t) = -2hfl F ( r , t) cos \urt - k .r + 

where F ( r , t) is the amplitude of the effective running laser beam, and 
where £Iq is the Rabi frequency of the atomic transition. 

An other important approximation is the rotating wave approximation (RWA), 
which consists in neglecting the off-resonant terms (i.e. with frequency u> + ll>q) 
in J3J). In our case, this approximation is supposed to be justified. Indeed, if 
the two effective levels are indirectly coupled through optical photons (as for 
Raman transitions for example), one can show that the Bloch-Siegert effect is 
negligible (4T1| . 

Finally, the evolution equation is equal to: 

- 1^ (t)) = iQ F [r op (t, to , tj ( e+l$o u p(Ml) Q J |^i (t)) 

where $ op is defined as: 

(f, h) = cut - uj (t - tO - "fc -Rop (*, <0 + 

The next unitary transformation corresponds to the usual passage into the 
rotating frame. In fact, there is an infinity of such transformations [I]. For 
example, the following family of transformations (indexed by the real parameter 
x): 

\<Pi{t))=U 2 (x,t,t 1 )\<p 2 (x,t)) (4) 

with: 

/ e -i4 oP (Mi)-a o 

U 2 (X,t,t 1 )=\ p -i<S> op (t,tl).(x-l) 



leads to the equation: 



Tt 1^2 (z, t)) = «M op (i, t) \ip 2 (x, t)) 



(5) 



with M op defined as: 



Afop (x, t) = 



x I u — ojq — k .R op — xS 



QqF 



r op 



(x — 1) w — luq — k .R op — [x — 1)8 



7 



where 8 (t,ti) is the "generalized recoil": 

5(t,h) = hkffk /2m 

and where the point above letters refers to the time derivative. 

Each value of x corresponds to a particular evolution equation. For example, 
the most "symmetric" choice is x = 1/2, whereas the most "physical" choice is 
x = 1, and leads to: 



M op (l,i) = 



A op i 



where A op i is the "generalized detuning" 0|: 



A op i (t, t\) — u) — cjo — k -R p — 8 

i.e., the operator which generalizes the usual "free" detuning in the presence 
of several gravitational, inertial and trapping potentials. It can also be written 
as: 



lip. 



Of) 



— I + aR op + f 



A op i (t, tx) = u> - ^ljq + k 

where P op is defined in the same way as R op jSJ 0] : 

P^ p (t, h) /m = C (t, t x ) ¥^ p + D (t, h) p^/m + ~$ (t, h) 



The expression of A op i (t, t\ ) can be easily interpreted considering the energy- 
momentum conservation for a non-excited atom absorbing a photon (u>, k ) at 
the instant t (the considered atom is then at the position R (t,to) with the 
momentum P {t, £q) when the absorption occurs): 



H ext (R(t,to),P (t,t ),tj+E b +hv = H ext (R(t,t ),P (t,t )+H~k,t)+E a 

which gives the non-operatorial version of the condition A op i(t,ii) = 
(exact resonance condition in the presence of the external potentials described 
by H ext ). 

This generalized detuning can also be expressed directly with the coefficients 
of Hext- For example, if H ext is constant, the first terms of its Taylor expansion 
(in a(t — ti) and j(t — ii) 2 ) are found to be 0]: 

A op i (t, t\) = u> — ujq — k £ — S — k .if (t — t\) — k .a.7\Z (6) 

m 

— 2k .a.—?- (t — ti) — k . (a 2 + 7) .t\1 (t — t\) — k .a.~g (t — t\) 2 — ... 
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where only the first five terms of the right hand side are non-negligible in 
usual experiments (weak rotations and acceleration gradients on the Earth). We 
can then use chirped laser pulses to eliminate the gravitational induced Doppler 
shift k . g (t — ti) and finally get back the usual "free" detuning. 

The other element of M op (x, t) which may depend on the two canonical 

operators is the effective amplitude F (^R op (t, ti) , tj . The two main sources 
of its spatial dependency are its transverse profile and the speckle due to the 
miscellaneous optical elements used to bring the lasers to the atoms. Generally, 
one can not neglect this speckle and the best is to map it, and then, to use the 
zones where the speckle is sufficiently weak. 

Concerning the laser transverse profile, one can show that it is seen roughly 
uniform by each individual atom of the initial atomic cloud (described by a 
statistical mixture of wave packets). One can therefore replace R op (t,ti) in 
F by its semi-classical action on a typical wave packet which evolves inside 
such beam splitters. For example, and as we see thereafter, R op {t,t\) can be 
approximated by: 

lC p (t, h) ~ A (t, h) ¥^ + B (t, ii) ( po + /m + ? (t, h) 

where and po are the initial central position and momentum of the con- 
sidered atomic wave packet. The main result is that F (^R op (t, ti) , tj ~ F (t, ti) 
is hence now independent of r^, and p^, (but is still time-dependent). 



4 Resolution methods 



The main problem in the integration of JSJ) is that, in the general case considered 
in this paper, M op (x, t) is a (2x2) matrix which depends both on time and on 
the two non-commuting canonical operators r op and p op . These are the two 
reasons why M op (x, t) does not commute with itself at different times, and why 
one can not apply the common rules to integrate JHJ) directly. 

However, in some particular cases, (JSJ) may depend on one time-independent 
operator only, and one can solve it analytically in the representation of this 
operator. It is thus important to list the maximum of these exactly solvable 
cases. This is the aim of the z(t) theory, initiated in j^I], improved in 02111] and 
generalized recently in [H] (for a detailed review, see [Hj and 0]). Among these 
exact solutions, let us underline the Landau-Zener model (solution with cylinder 
parabolic functions) which accounts for the effect of a time-independent and 
uniform acceleration during the atomic splitting, and the Rosen-Zener model 
(solution with hyperbolic secant functions) which is significant in the study of 
matter wave solitons. 
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Some other analytical methods are particularly useful to deal with the equa- 
tion J5j): Floquet theory 03| f° r periodic time-dependence and its generaliza- 
tions (multi-periodic Floquet method (t,f) theory 05]...); bands theory 
(i.e. use of Bloch states) when one can not make the RWA 015]; use of quasi- 
probabilities (Wigner and Shirley representations) and phase space functions 
(Q function of Hushimi and Kano, P distribution of Glauber and Sudarshan) 
when some QED effects can not be neglected (for a recent review, see 03)- ■ ■ 

Apart from these particular exactly solvable cases, it is always possible to 
write the general solution of © as a formal development. This development 
may be linear (Dyson) or not (Magnus, Fer, Cayley. . . ), and may preserve the 
unitarity (for a review of the recent advances concerning the Magnus expansion, 
see 0U0S1E])- However, due to the entanglement of operators and p^, in 
the different terms of these developments, it is impossible to choose any repre- 
sentation which leads to an analytical expression of the solution. This problem 
can be solved either in eliminating one of the two canonical operators directly 
in the equation (JSJ) ( "operatorial elimination method"), or in approximating the 
generalized detuning A op i (t, ii), or finally in solving the equation numerically. 

The operatorial elimination method, which is detailed in 4 , leads to a double 
development, easily calculable but rather long, which is why its use would be 
limited to the numerical domain. As for the approximations, one has already 
underlined that the effect of rotations and acceleration gradients may often be 
neglected in the generalized detuning (leading to a trivial integration of ©). If 
not, several tactics may be used 0]: 

1. freezing A op i (t, t±) at a particular time (mid-time for example) or taking 
a temporal average, and choosing the resulting time-independent repre- 
sentation; 

2. or replacing and p^, by their semi-classical value (WKB approxima- 
tion); 

3. or replacing only one of these operators (r op for example) by its action 
on the initial atomic wave packet (or on a wave packet closer to the final 
solution) 

When the equation (0 is made scalar, one can use either the previous an- 
alytical methods, or the previous developments (and truncate them when they 
converge, see [50] 021), or some intermediate approaches which are based on the 
eigenstates of the matrix M op (x, t) (super-adiabatic scheme jSJ or successive 
adiabatic states method 0]). The latter are extremely interesting because they 
directly provide important information on the solution, like the true energies of 
the system "atom - laser - other external fields" and the corresponding group 
velocities. 

Several numerical methods may also be implemented (for a recent review, 
see |52p: Magnus expansion, median exponential rule, Runge-Kutta method, 
Strang-Marchuk- Trotter method, Chebyshev or Lanczos approximations. . . 
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Finally, we obtain a solution which is more or less close to the exact solution 
of ©: 



\fsol (*)) = S op (t, t ) |<^2 (*o)) - |^2 (t)> 

with the evolution operator: 

S op (t, t ) ~ T ^exp ^ jf M op (t') dt' 

which can be written explicitely as a S matrix between the initial and final 
atomic states: 

f \b S ol(t)) \ _ f Sbb,op{t,to) Sba,op{t,to) \ ( |^2 (to)) \ 
V \ a sol (t)) J \ Sab, op (t, to) Saa,op(t,t ) J \ \a 2 (to)) / 

The resulting action of this S matrix on the initial atomic wave packet is a 
possible change of internal state (generalized Rabi oscillations) and a structuring 
into several wave packets that can be quite different from the initial one. In 
particular, the group velocities of the created wave packets may be identical 
or not (Borrmann effect, see below) according to the value of the generalized 
detuning. 



) 



5 Return to the initial picture and equivalent 
instantaneous interaction 

5.1 Solution in the initial representation 

Once we have obtained the matrix 5 , op (t,to), we can do the reverse unitary 
transformations of and (J2J, and go back to the initial representation. One 
obtains (written here for the symmetric choice x = 1/2): 

\^ sol (t)) = e- iI *o 5 ( t '' t ^ dt ' / \U 1 (t,t 1 )U 2 (l/2,t,t 1 ) (7) 
So P (t, to) U^ 1 (1/2, t , h) [7f 1 (t , h) |* (t )) 

which expresses the link between the general solution \^ so i (t)) and the initial 
bet |* (to)). 

It may be noticed that the solution depends on tx, the arbitrary time in- 
troduced to define the previous unitary transformations. This instant has no 
physical meaning and can be removed explicitly at each step of our calculations 
J3|. However, it is more interesting to keep it and eventually assign it a particu- 
lar value (the central time of the interaction for example) which may be useful 
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both for calculations and the interpretation of the obtained solutions. The rel- 
evance of this instant to describe the atomic beam splitting as an equivalent 
instantaneous interaction (generalized ttt scheme) will appear more clearly at 
the end of this part. 

It is easy to interpret the expression of \Q S oi (fit) by writing the typical terms 
exp (±|$ op (t,to)) of U^ 1 as a product of exponentials. Indeed, let us consider 
the following initial ket: 

'*«»» = ( M<o» 

which represents the incident atomic wave packet (atoms in the lower state 
a). It results that the upper component of \^ S oi (t)) is: 

e l8lba U ext (t,h) e l~k A (t,ti)^p e fkB{tM)p?p/™ (g) 
S ba ,o P (t,t ) e f^(*o,ti)^ e itB(t , tl )p^/m Uext (tl,to) \a(t )) 

which can be interpreted as follows: 

1. evolution from t to ti due to U ext , i.e. to the gravito-inertial and trapping 
potentials described by H ext 

2. translation of —B (to,t\) |— - in position and +A (t , t\) 4p in momentum 
due to the term exp ^| k A (to, t\)r^j . exp ^| k B (to, ti)p^,/rnj 

3. action of the non-diagonal element of the S matrix Sb a ,op (t,to) (change 
of internal state and dispersive structuring) 

4. again two translations: — B (t, t\) in position and +A (t, t\) in mo- 
mentum 

5. then the external evolution from t\ to t due to H ext 

6. and finally a phase factor exp (i6n, a ) with 



E b (t - t ) ( w ~ (* - *o) - (T (t, h) + T {to, h)j 

9lba = 7 Uto — 



+ [A (t, h) B (t, h) + A (t , h) B (to, ti)) ~k-f 
Similarly, the lower component of l^soi (t)) is equal to: 



S J^lldt> 

to 4 



e 



iOi, 



a U ext (t,t\) k A{tM)r^i e -% k B(tM)p^t/m j-g-j 



with the following similar interpretation: 
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1. evolution from to to t\ due to H ext 

2. translations of — B (to, ii) §^ in position and of +A (to,t\) in momen- 
tum 

3. action of the diagonal element S aa . op (t, to) of the S matrix (dispersive 
structuring without internal change) 

4. translations of +B (t, t\) 4^- in position and of —A (t, t\) 4p in momen- 
tum 

5. evolution from t\ to t due to H ext 

6. and a phase factor exp (i9i aa ) with 



E a (i - to) ( w ~ ^ ^ - to ^ ~ k ( ^ ^ ~ ^ 4l V 



'laa 



h 2 



which differs from #ib Q by urf + <p — k . £ (£, ti) . 
5.2 Generalized atomic Borrmann effect 



In certain conditions, the transfer matrix S op may have only a weak effect on 
the external state of the incident atoms. In this case, the center of the two main 
wave packets (associated with the upper and lower adiabatic atomic states, see 
below) evolves along the same trajectory during the splitting: it is the (atomic) 
Borrmann effect. 

This effect is well known in dynamical diffraction of X rays |53| , of neutron 
waves and more recently of atomic waves ^U]. Historically, this effect 

was defined as a more specific phenomenon, valid for any kind of waves diffract- 
ing in an absorbing crystal. It can be stated as: "for a certain angle of incidence 
with respect to the crystal surface, the propagation of the incident wave packet 
inside the crystal is made with no attenuation along a unique trajectory which 
is orthogonal to the crystal surface" . This particular angle of incidence is the 
well known Bragg angle defined as: 

P. (po" + P/2) = 

where po is the central momentum of the incident wave packet and P is the 
quantum of momentum which is communicated to the diffracted wave packet 
(the Bragg condition is written here for the first order diffraction). 
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Conversely, two wave packets with two different trajectories are created in- 
side the crystal if this condition is not fulfilled (defining the well known "Bor- 
rmann fan"). Furthermore, this anomalous transmission is very sensitive to the 
Bragg condition, and any deviation from it greatly amplifies the angle between 
the two trajectories [To] . 

In the case of atom-laser interactions, there is no absorption of atoms and 
the Bragg condition (which means nothing but the energy-momentum conser- 
vation, i.e. the resonance condition) is partly relaxed due to the atomic internal 
structure: 

lu — luq — k . ^p~Q + hk /2j jm = 

To date, this atomic Borrmann effect was only studied in the free case [TP] 
(for which H ext is limited to the usual kinetic Hamiltonian p 2 /2m) or in the 
presence of a time-independent and uniform acceleration 12 . In this paper, 
it is obtained in the presence of the various gravitational, inertial and trap- 
ping potentials described by H ext . By the way, we will show that the common 
Borrmann trajectory is equal to the average of the two extreme trajectories: 
A ro + B po/m + £ (atom absorbing a photon at the final time t) and A 
¥q + B ( po + h k^j /m + £ (atom absorbing a photon at the initial time to)- 

Indeed, let us consider the previous example (atoms initially in state a) and 
suppose that S op has a negligible effect on the central position and momentum 
of the corresponding incident wave packet. Then, according to the expression 
(J7J and the Ehrcnfcst theorem, we can show that the central position ro of the 
initial wave packet is changed into: 



f*c (t, h) = A fr^- B^-j +B\p5+(l + A) /m + ? 
where po is its initial central momentum. 

Finally, thanks to some simple properties of ABCD matrices, we obtain the 
previously stated result: 



rc (t, t ) = ATo* + B [ $ + ) /m + f 



2 



Ar\l + Bp^o/m + T) + (Atq + B (po + hkj fm + ( 



This unique central trajectory differs by B (t, to) -H k /2 from the one ob- 
tained without any splitting potential. It means that, even for the atoms which 
are finally in the same internal state as the initial one (no effective internal 
change), there is a non-trivial change of their central trajectory, which results 
in a measurable spatial shift at the end of the interaction (see Figure^. 
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Outgoing 




Figure 1: The continuous lines represent the central atomic trajectories (group 
velocities); the dashed lines refer to the equivalent instantaneous interaction (ttt 
scheme); and the dotted line is the extension of the initial atomic trajectory 
(infinitely thin modeling). 




As for the central momentum, we obtain similarly the two following mo- 
menta: 



p* c (t, t ) = mCn, + D 



which differ from each other on h k , as expected. 

On the contrary, S op may have a non-negligible effect on the external state 
of the incident wave packet (change of central position and momentum, respec- 
tively due to the addition of a non-trivial group velocity and to the change of 
momentum distribution). In this case, one can show that the initial wave packet 
is split into two main wave packets, which evolve along two different trajectories 
(i.e. with two distinct group velocities) which form the atomic Borrmann fan. 
For each of these wave packets, we can do the same calculation as before and 
obtain the trajectories in the initial frame 0]. 
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5.3 Generalized ttt scheme 



However, in both previous cases, it is noticeable that the expression (JJJ can be 
written as a product of three evolution operators: 

l^soi (*)) = Ui (t,t 1 )S 1 (r^,p^,,t,tQ,ti)Ui (h,t ) |* (*o)> 

where U\ (ti, to) and U\ (t, ti) describe the evolution due to H + H ext only, 
and where S x represents the evolution part depending on the splitting poten- 
tial V. The aim of this arrangement is to clearly separate the effect of H ext 
and V, and to describe the interaction as an effective instantaneous interaction 
(generalization of the ttt scheme introduced in |14p. 

In addition to the link with the infinitely thin modeling of atomic beam 
splitters, we aim at providing a clear and practical beam splitter modeling, in 
the presence of the external fields described by H ext , which can be used easily 
in atom interferometric phase shift calculations 

To clearly separate the splitting terms from the translation and phase terms 
in (jSJ and we can transform the following expression: 



»(*o,ti) c =r-*\ „Ti*o P (*o,ti) 

U UV \' Op 7 fop J ° 



S U v r^T B(t ,ti) —,po^±A(t Q ,ti 



h k v ~ , . > hk 

2m 



thanks to the following algebraic properties: 
e A .f(B).e- A = f(e A .B.e- A ) 

e A .B.e- A = B + [A, B] + i [A, [A, B]] + ... 

where A and B refer to operators or square matrices, and where / is a function. 
Finally, the diffusion matrix Si can be written as: 



Q _ ( ^1,66 Si, ba \ 
I Si, a ft Si, aa J 



where its elements are equal to (a and b are the lower and upper states 
respectively) : 

1. for the a — ► a transition: 

Si,oa = e^ei^-^e-^-^Saa (r£-B (* ,*i) ^,p£ +A(t , h) - 
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with: 



A(t,h) -A{t a M) 



hk 



B (t,h) -B(t ,h)hk 



m 



<I> = 



+ 



1 

+ 2 

hk 
8^" 



(w - w ) (t - t Q ) - ~k. (j (t, h) - ~t (to, *l)) 
i4 (t, ti) B (t, ti) + A (t , ti) 5 (to, ti) - (t , ti) B (t, ti)| fc* - / 
2. for the a — ► b transition: 



hk 



op - B (to, ti) ~^>P°p + ^ (*0)*i) 



with: 

= + ^(Mi)+/(*o,*i)^ 



B(t,ti)+B(to,ti)nfc 



t + t 

w — CJQ 



m 
t + t 



ti + <A - fc . 



£ (t,ti)+ £ (t ,ti) 



8^ 



A (t, ti) 5 (t, ti) + A (t , h) B (t , h) + 2A (t , h) B (t, h) 
3. for the b — ► a transition: 

Si,„ 6 = e^ei^-^e-i^-^Sab (t^ p + B (t ,h) ^,P? P - A {to, h) 



t6 -^ldt' 



to 



hk 



with: 

p2 = _*(Mi) + A(«d,ti)^ 



r a 6 = + 
= + 



B(t,ti) + B(toM)hk 



m 



t + to 

OJ — LUq 



e (t,ti)+ e (t ,ti) 



+ 



hk 
8m 



A (t, ti) B (t, h) + A (t , ti) B (to, ti) + 2A (t , ti) B (t, ti)l ~k - f 

J t 



^(t'M) dt , 



17 



4. for the b — ► b transition: 

Si, 66 - e^efrZ-rZe-tt'^Su (r£ + B (t , h) ^,P? P ~A(t , h) ^ 
ou: 

Pbb = H ^ " fc 

B(t,ti) -5(f ,ti) 



»'6fe 



l r 

~2 

8m 



(w - uq) (t - t Q ) - ~k. (? (t, ti) - 7 (*o, ti)) 
A (t, ti) B (i, ti) + A {t a , h) B (t , h) - 2A (to, ti) B (t, h) 



The interpretation of these terms is simple and constitutes the core of the 
generalized ttt scheme: 

1. from to to t\, the initial ket \*f? (to)) evolves with Hq + H ext only (as if 
there was no splitting potential); 

2. at t\, this evolved ket is subject to an effective instantaneous interaction 
which modifies its structure (structuring into several wave packets, see part 
and the illustration of the ttt scheme on figure [2J , its central position 
and momentum, and its phase; 

3. from ti to t, the obtained wave packets evolve with Ho + H ext only, as 
before t\. 

Eventually, we can show that the t\ value which simplifies the most the 
previous expressions is the central time of interaction: 



f i = t 



1/2 



t + t 



2 

Let us see now what are the features of the wave packets which appear inside 
the beam splitters. Before studying this wave packet structuring in the general 
case (i.e. in the presence of trapping and gravito-inertial potentials), we will 
focus on the simple "free" case (i.e. without any such other external potential). 



6 Structure of the solutions in the free case 



The study of the free case is important for at least two reasons. First, as we have 
seen in partOH the effect of the quadratic terms of H ext (rotations, gradients, 
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Outgoing 




2 



Figure 2: Atomic beam splitting in the free case and for an inelasticity parameter 
different from zero. The dashed lines refer to the ttt scheme, whereas the dotted 
line represents the common atomic trajectory of the Borrmann effect. 



trap. . . ) can often be neglected during the resolution of equation (J5J, and any 
uniform acceleration can be easily removed by modulating the laser frequency. 
In this case, the equation JSJ is equivalent to the one stemming from the free 
case. Second, and as far as we know, to date, there has been no global study of 
matter wave beam splitters even in the simplest free case. Indeed, each of the 
previously quoted works aims at studying only a few aspects of the atomic beam 
splitting: mechanical effect, internal and external splitting, group velocities. . . 
Actually, the question consisting in finding what exactly goes out of an atomic 
beam splitter is still relevant, even for a simple laser potential (temporal square 
amplitude) and an incident Gaussian matter wave packet. 

It is therefore necessary to specify what is the exact structuring of this 
incident wave packet, namely the number of created wave packets (two main 
ones for each transition amplitude), their amplitude (Rabi oscillations, velocity 
selection, anomalous dispersion. . . ), their group velocity (Borrmann fan), their 
central momentum and their phase, and eventually how these quantities evolve 
in time. 
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6.1 Group velocities and adiabatic states 



The laser-atom interaction induces particular states, the adiabatic or dressed 
states, which are nothing but the eigenstates of the interaction. If F (t, t\) is a 
unit amplitude square pulse between to and t, these adiabatic states are simply 
the eigenvectors of the matrix iM (1/2, t). Therefore, the two corresponding 
eigenenergies are: 



Te(7) = Th^nl + [A l (]?)/2r 

with A; ( ~p) — lo — ujq — k . ~p/m, and the corresponding group velocities 
are simply obtained by deriving these dispersion relations with respect to the 
momentum: 

±Vg(P)=T- 



2 2m 



where y (~p) is the well known "off Braggness parameter" introduced in 
neutron optics [SlIKj and defined here as: 

V (~P ) = (u — uio— k .Jj/rnj /2fL 

As we can see on (JHJ and 1(5} . the matrices S U v,op act on a wave packet which 
is shifted by a global momentum of +H k /2 or —h k /2, depending on the initial 
internal state of atoms. For the example previously described (atoms initially 
in the lower internal state), we obtain the following parameter: 

y{p+ hk/2) = y+ (~p) = (w-wo- ~k-~p/m - S^j /2f> (10) 

which can be called "inelasticity parameter" as it refers to the way the 
resonance condition is fulfilled In the initial representation, these group 
velocities become: 

+ — (i± y+(][) 

m 2m 1 r. ' > x 2 



\Jl + y+ (~pY 



The examination of these velocities leads to several important conclusions. 

First, the difference between momentum and group velocity in the presence 
of an electromagnetic field is naturally confirmed. For a weak inelasticity pa- 
rameter \y + (p)\ <C 1, we obtain only one group velocity for both the adiabatic 
states (atomic Borrmann effect, see part 15.2(1 : 

~p hk 
m 2m 
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whereas for \y + (~p)\ 3> 1, we obtain the two extreme velocities (defining the 
Borrmann fan): 



P , P 
— and — 



hk 



rn 



But in this case, the beam splitter is inefficient as we will see thereafter. 

However, we can show |lj that these group velocities are closely linked to the 
average momenta of the considered two level system, and that they are more 
precisely equal to the most probable momenta of adiabatic states (divided by 
m). 

For y 7^ 0, two distinct atomic wave packets are created in the beam splitter, 
and their physical separation may be observable in certain conditions (more than 
few /itm after 10 -4 s for an initial atomic coherent state of 1 /im width). 

Finally, it is noteworthy that these group velocities depend on p . Therefore, 
the (optical) medium where atoms evolve is dispersive and this effect leads to 
the phenomenon of anomalous dispersion (see part I6.4|l . 



6.2 Rabi oscillations 

Let us consider the free solution of equation JSJ . If the internal state of initial 
atoms is the lower state, we obtain in the initial picture the following lower 
state: 



with: 
T aa (~P,t) = cxp 



T aa (P^>,t) \a(t )) 



+ — /2mh+- —2.\( t -to) 




and the following upper state: 

— il ujt — k 



T ba (p^>,t) \a(to)) 



with: 



Tba ("p , t) = exp 



, Eh I _» hk \ . . 8 uj — cjn i 
'l + r + T /2mh+- + ^-^\(t-t ) 




S aa and Sba having the following usual expressions: 



S aa (~p, t-to) = COS [e ("p ) . (t - to)] - % 



y{~p. 



:sin [e (p) . (t - t )] 
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Sba (p, t — t a ) = i , = sin [e (p) . (t - t )] 

\]^ + y(~p? 

The temporal evolution of the initial state is therefore sinusoidal with an am- 
plitude equal to 1/ (l + y + (~p ) 2 ^j and a frequency equal to flo\J 1 + y+ ( ~p) 2 /2tt. 

Tuning the interaction parameters, we can realize true atomic mirrors or "7r 
pulses" (transfer of all the atoms from one state to another) and semi-reflecting 
plates or "7r/2 pulses" (50-50 splitting), whose efficiency depends in a crucial 
way of the value of y+ (~p) which is actually taken by the central momentum of 
the incident wave packet. 



6.3 Velocity selection 



This effect comes from the fact that the two pre-factors 1/y 1 + y + (~p*) and 

V+ (~p) /\A + V+ Vp ) 2 - m the expression of S aa and Sb a , depend on the mo- 
mentum p , and more particularly on its part which is collinear to the laser wave 
vector k (transverse momentum). Therefore, the terms S uv act as momentum 
filters on the incident atomic momentum distribution. 

For example, in the case of the a — > b transition, one obtains the well known 
sinus cardinal filter. It is characterized by a central lobe with an amplitude of 

sin {Qoi~} and a total width of 4£!o {it /VLqt) 2 — l^j jk (for a fixed interaction 

duration of r = t — to). This width may be less than the transverse width of 
the incident atomic momentum distribution. The resulting (transverse) velocity 
selection is very useful in atom interferometry (to increase the fringes contrast) 
and constitutes the basis of Raman cooling |57| . 

Furthermore, this central lobe is centered on the transverse momentum pf. 

pi = m (lo — ujq — 6) jk 

which can be quite different from the initial transverse central momentum 
Pq. k jk. If so, the central momentum p^ of the filtered distribution is distinct 
from p$ and is bounded by po and pi . Consequently, the central momentum of 
the outgoing wave packets (in excited state) is not p^ + H k , as one would think 
in view of the mechanical recoil effect of light, but pi + hk (see Figure [3J|. 

Apart from the central lobe, the sine filter has nodes and sidelobes with 
an amplitude that rapidly decreases (this particular structure is linked to the 
pulse shape, and can be softened or removed by the use of apodization functions 
[Blackman for example] to tailor laser pulses). If the incident atomic distribution 
is sufficiently broad to encompass one or several sidelobes, the filtered matter 
wave packet will then be structured into several wave packets whith central 
momenta quite different from po, pi or p c . 



22 




With respect to the other transition amplitude (a — ► a, i.e. without change 
of internal state), we obtain the complement filter and the structuring can be 
studied in the same way 0]. 



6.4 Anomalous dispersion 



This effect refers to the modification of spreading of atomic wave packets inside 
a beam splitter. It is more convenient to examine it for an incident atomic wave 
packet which has a thin momentum distribution, although it can be studied and 
modeled for any momentum distribution (see 0]). 

As we consider the wave packets which evolve inside the beam splitter, we 
have to go into the adiabatic states. In the adiabatic picture, these wave packets 
have the two energies (written here for the a — > b transition): 

±m y/l + y+ (^) 2 

where the term \J \ + y + {~p) 2 can be Taylor expanded with respect to 
k . [~p — po) /2t7iS1o as: 
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T 2mO (l + yl + f 2 4m«lg 
where yo+ is defined as: 

2/o+ = 2/+ (po) 

and where 8 is the complement matrix of the recoil term 5: 
=> hk/k 



2m 




In the initial picture, the first order term of this expansion gives the group 
velocities obtained in part 16.11 As for the second order term, it corresponds 
to an additional dispersion and it indicates that one wave packet spreads more 
than the natural spreading, and that the other one spreads less. In certain 
conditions, this spreading can even be stopped or changed into a contraction 

In the case of a non-thin incident atomic wave packet, the main results of this 
study are still valid provided yo+ is changed into y c + = y+ (p~c)- One can then 
model the outgoing atomic wave packets. A simple but powerful way to do this 
is the Gaussian modeling which consists in writing these wave packets as Gaus- 
sians. This "strong field ttt modeling" , and more generally the use of Gaussian 
wave packets, is found to be particularly relevant in atom interferometry (see 

m 



7 Structure of the solutions in the general case 

We have already seen in part 0] how to deal with the double non-commutation 
problem which appears in A op i and in the resolution of In particular, we 
have seen that one can often neglect, in the expression of A op i, the terms de- 
pending on the quadratic terms of H ex t, namely a, j3 — 1 and 7. In certain 
configurations (when k is orthogonal to g or when lu is modulated to compen- 
sate the gravity induced Doppler shift), the generalized detuning is even equal 
to the free one and it boils down to the free case in (and only in) the resolution 
of 0. 

In this case, we obtain the same previous adiabatic energies and consequently 
the same adiabatic group velocities, Rabi oscillations, velocity selection and 
anomalous dispersion effect as stated before. In the initial picture (i.e. in the 
lab frame), the elements of the previous Sj matrix are therefore expressed as 
(u, v — a,b): 



24 



Si,„„ = exp [i$uv] ■ exp 





i > > 




i 






exp 


^7 Top • Puv 


. exp 


rn Pc 


p-Tuv 


■iJuv 1 Pop ' O / 



where the elements S uv are equal to the ones obtained in the free case (see 
the previous part). 

Finally, the only differences from the free case lie in the expression of: 

1. the effective incident wave packet at time t\ (which is equal to the initial 
wave packet evolved from to to t\ thanks to Hq + H ext ). In particular, 
its central position and momentum are no more ro + Po (*i — to) /m and 
Po but R (ti,to,ro,Po) and P (ti,to> ^o>Po) (expressed with the ABCD 
matrices, see part|2J) 

2. the terms r\^, pJa and $ m which are detailed in part 15.31 

3. U\ (t, t\) which accounts for the evolution from t\ to t due to Hq + H ext . 

Then, we can do the same Gaussian modeling as in the free case previously 
studied (see 0j). 

In some cases however, we can not neglect the effect of (unavoidable) grav- 
itational, incrtial and trapping potentials for the resolution of JSJ), and it is 
important to examine what are the changes, due to H ext , of the properties 
listed before (group velocities, Rabi oscillations, velocity selection and anoma- 
lous dispersion). Moreover, we have already underlined that the key parameter 
is the inelasticity parameter y which is proportional to the generalized detuning 
A. This y parameter generally depends on the two canonical operators and 
Po^,, but in the case of linear potentials (uniform acceleration for example) it 
depends on one operator (p op ) only. 

Before dealing with the general case, let us consider the effect of a uniform 
(but time-dependent) acceleration g (t). From the equation , which becomes 
scalar in momentum representation, one can extract the adiabatic energies and 
then the adiabatic group velocities (written here for t\ = to). 



±Vg~0 (*,*o) = T" 



yo hk 



with (F is taken constant and equal to 1 for simplicity): 

rt 



yo = y(po) = yofr 



g (t') dt , /2tt 



09) 



where yo free is the inelasticity parameter obtained without g (free case, see 



In the initial picture, this result becomes eventually: 
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m 2m J tQ 

where two sources of atomic trajectories bending can be identified: the com- 
mon gravitational bending coming from the third term J t g (t')dt', and an 
"anomalous gravity induced Doppler bending" coming from the time-dependent 
adiabatic group velocities. As for the free case, this anomalous bending can be 
described either in terms of effective mass tensors ^3 ^2 or m terms of position 
(and time) dependent effective refractive index inside the beam splitter ^2 [55] . 

Let us remark that the acceleration induced Doppler term — k .~g may pro- 
duce a non-trivial bending of central atomic trajectories. Indeed, if yofree and 
the scalar k . g are non zero, and if yofree has the same sign as k . g , then the 
trajectories bending due to the k . g term makes these trajectories get closer 
(as usually considered in literature). On the contrary, if yo/ree has the opposite 
sign of k . g , then these two main trajectories repel each other (see Figure 0}. 
In both cases, if the vectors k and ~g point at the same direction (k .~g > 0), 
downward for example, then one of the two main atomic wave packets will al- 
ways be accelerated upward, even if all "forces" (acceleration and laser pulse) 
seem to push the atoms downward. 




Figure 4: Central atomic trajectories inside a (temporal) beam splitter in 
the presence of Earth gravity: adiabatic picture. Numerical application: 
k = 10 7 to-\ g = 10m.s~ 2 , n Q = 3.10 3 s-\ y ofree = ±0.1, t = 0s. The 
central atomic trajectories either attract or repel each other. 
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Furthermore, this particular behavior is not limited to the adiabatic picture 
and can be observed fully in the laboratory frame (i.e. in the initial picture). 

/ \ 3/2 

Indeed, if d > 2fl ( 1 + yof ree ) (which corresponds, for yof re e = 0, to the 

"Bragg regime" , as it is defined in in contrast with the "channeling regime" 
for which 5 < 2Qq), this anomalous upward acceleration exceeds the downward 
acceleration of gravity (for t — to much less than l/f^o)- As a result, some atoms 
are accelerated upward in the lab frame, even if the action of g and laser 
beams is directed downward (see Figure |3J. Of course, this anomalous "anti- 
gravitational" bending can be explained simply by considering the conservation 
of momentum and energy during the interaction process. 




Figure 5: Central atomic trajectories inside a (temporal) beam splitter in the 
presence of Earth gravity: initial picture (lab frame). Numerical application: 
k = lO 7 ™- 1 , g = 10m. s- 2 , n„ = 3.10 3 s _1 , y free = -0.1, t = 0s. For short 
times, the anomalous upward acceleration exceeds the downward acceleration 
of gravity. 



In the general case of an at-most-quadratic H ext (with quadratic terms like 
rotations, gradients of acceleration, trapping potentials, etc), the generalized 
detuning A op i depends on the two non-commuting canonical operators, and it 
can not be made scalar in any representation. As we saw previously, one of 
the most relevant strategies is to approximate, in the resolution of equation (J5J , 
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the effect of one of these two operators by considering its action on a typical 
atomic wave packet which evolves inside the beam splitter. In the free case, 
the Gaussian approximation of these typical wave packets, called here "wp" for 
convenience, leads to (in the adiabatic picture): 



r 



op 



(wp) ~(r$±£ (f, t ) dt' + (U T U AD ) ■ (~P ~ Po)^ .wp (11) 



where Uq is linked to the complex momentum width of the initial atomic wave 
packet, and where Uad accounts for the anomalous dispersion phenomenon. 

In fact, we can make Uad arbitrarily small (for a sufficiently small interaction 
duration) and incorporate the time-independent matrix Un into the expression 
of y (see © and Qll)|0. and consider only the two first terms of l|llfl . Finally, 
the quantity ~x (t) = f (t' , t ) dt' meets the following first order differential 
equation (which can be solved numerically): 

d_ yt-~k.A(t ) t ).l?/2n F(t) jf^r.hk 

dt , 71 y u; 2m 




k .A(t,t ) .~x/2tt F{t) 



with: 



VE= \u-ua-6- k.A(t,to).ro - k.B(t,to).po/m- k.£ (Mo) /2QqF (t) 



A simpler approximation (WKB approximation) amounts to neglecting also 
It ^9° (^'' *°) in the expression fTTft . In this case, we obtain readily: 

y , . Tin — hk 

±^ (i,f ) = T _^_B(i,i )_ 

which gives a good approximation of the group velocities inside a matter 
wave beam splitter when this latter is subject to the various external potentials 
described by H ext . 

For example, for a non-uniform (but time-independent) acceleration (3,7) 
(or in the case of a trapping potential, for which the sign of 7 has to be reversed), 
one obtains: 



2/o = 

- k 



(uj — ujq — S — k . cosh (7 1 / 2 (t — to) J ■Pa/ m 
7 1 / 2 sinh (7 1 / 2 (t - t )) .rj - ~fc . 7 " 1/2 sinh (7 1 / 2 (t - t )) .~g) /^ F (t) 

Similarly, in the case of a (time-independent) rotation f2 , one obtains: 
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yo 



(lj — uj - 5 - k .1Z (t, t ) . po/m — Q x (ro +pa (t — to) /m) ^ /2Q Q F (t) 



where the rotation matrix lZ(t,to) can be seen as acting either on the atomic 
wave packet or on the wave vector k . 



8 Conclusion 



In conclusion, we have shown in this paper how to solve the problem of matter 
wave splitting in the presence of various gravitational, inertial and trapping 
potentials. In particular, we have seen how the resonance condition between 
the splitting potential and the effective two-level atoms has to be changed. 
Then, we have shown how to express this triple interaction "matter - splitting 
potential - other external potentials" as an equivalent instantaneous interaction 
(generalized ttt scheme). 

Finally, we have investigated in detail what is the dispersive structuring of 
an incident atomic wave packet inside such beam splitters, both in the free 
case (for which H ext is reduced to J? 2 /2m) and in the general case of H ext . 
Several significant features of the solutions have been studied: group velocities, 
generalized Rabi oscillations, velocity selection, anomalous dispersion effects... 
In the light of this study, the generalized ttt scheme leads to a very practical 
and efficient (Gaussian) modeling of atomic beam splitters which is particularly 
relevant for atom interferometric signal calculations E] ■ 

It is worth pointing out that these results stem from a strong field theory 
for both the splitting potential and the other external potentials described by 

Hext- 

However, several points still have to be cleared up: for instance, the prob- 
lem of the effective mass change which occurs when the atomic internal state is 
changed, and which leads to non-trivial (small) relativistic corrections. Further- 
more, it is necessary to extend our formalism to additional external potentials 
which are more than quadratic (in position and momentum) if we want to in- 
vestigate the effect of van-der-Waals, Casimir or Yukawa-type potentials on the 
matter wave splitting. More generally speaking, it would be interesting to go 
beyond the various approximations listed in part [21 and in particular beyond 
the two-beam approximation. 



References 

[1] T.L. Gustavson, A. Landragin and M.A. Kasevich., Class. Quantum Grav. 
17, 2385 (2000); A. Peters, K.Y. Chung and S. Chu, Metrologia 38, 25 



29 



(2001); J.M. McGuirk, G.T. Foster, J.B. Fixler, M.J. Snadden, and M.A. 
Kasevich, Phys. Rev. A 65, 033608 (2002); G. Wilpers, T. Binnewies, C. 
Degenhardt, U. Sterr, J. Helmcke, and F. Riehle, Phys. Rev. Lett. 89, 
230801 (2002). 

[2] Ch. J. Borde, N. Courtier, F. du Burck, A.N. Goncharov and M. Gorlicki, 
Phys. Lett. A 188, 187 (1994). 

[3] Ch. Antoine and Ch.J. Borde, J. Opt. B: Quantum Semiclass. Opt. 5, S199 
(2003). 

[4] Ch. Antoine, Contribution a la theorie des interfcrometres atomiques, 
Ph.D. thesis (in French), Universite Pierre et Marie Curie, Paris (2004). 

[5] F.T. Hioc and C.E. Caroll, Phys. Rev. A 32, 1541 (1985). 

[6] K.-A. Suominen and B.M. Garraway, Phys. Rev. A 45, 374 (1992). 

[7] J. Ishikawa, F. Riehle, J. Helmcke and Ch.J. Borde, Phys. Rev. A 49, 4794 
(1994). 

[8] L. Carmcl and A. Mann, Phys. Rev. A 61, 052113 (2000). 

[9] A.M. Ishkhanyan, Opt. Comm. 176, 155 (2000). 

[10] M.K. Obcrthalcr, R. Abfaltcrcr, S. Bernet, J. Schmicdmaycr, and A. 
Zeilinger , Phys. Rev. Lett. 77, 4980 (1996); Ch.J. Borde and C. 
Lammerzahl, Ann. Phys. (Leipzig) 8, 83 (1999). 

[11] B. Eiermann, P. Treutlein, Th. Anker, M. Albiez, M. Taglieber, K.-P. Mar- 
zlin and M. K. Obcrthalcr, Phys. Rev. Lett. 91, 060402 (2003); B. Eier- 
mann, Th. Anker, M. Albiez, M. Taglieber, P. Treutlein, K.-P. Marzlin and 
M. K. Obcrthalcr, Phys. Rev. Lett. 92, 230401 (2004). 

[12] C. Lammerzahl and Ch.J. Borde, Gen. Rcl. Grav. 31, 635 (1999). 

[13] C. Lammerzahl and Ch.J. Borde, Phys. Lett. A 203, 59 (1995); K.-P. Mar- 
zlin and J. Audrctsch, Phys. Rev. A 53, 1004 (1996). 

[14] Ch.J. Borde, Gen. Rcl. Grav. 36, 475 (2004). 

[15] M.A . Home, Physica A 137, 260 (1986). 

[16] H. Rauch and S.A. Werner, Neutron Interferometry (Clarendon Press, Ox- 
ford, 2000). 

[17] O. Carnal and J. Mlynck, Phys. Rev. Lett. 66, 2689 (1991); F. Shimizu, K. 
Shimizu and H. Takuma, Phys. Rev. A 46, R17 (1992). 

[18] D.W. Keith, M.L. Schattenburg, H.I. Smith, and D.E. Pritchard , Phys. 
Rev. Lett. 61, 1580 (1988). 



30 



[19] J.F. Clauser, Physica B 151, 262 (1988). 

[20] F. Shimizu, Phys. Rev. Lett. 86, 987 (2001); T.A. Pasquini, Y. Shin, C. 
Sanner, M. Saba, A. Schirotzek, D.E. Pritchard, and W. Ketterle, Phys. 
Rev. Lett. 93, 223201 (2004). 

[21] Yu.L. Sokolov, Sov. Phys. JETP 36, 243 (1973). 

[22] J. Robert, Ch. Miniatura, S. Le Boiteux, J. Reinhardt, V. Bocvarski, J. 
Baudon, J. Europhys. Lett., 16, 29 (1991); Ch. Miniatura, J. Robert, S. Le 
Boiteux, J. Reinhardt and J. Baudon, Appl. Phys. B 54, 347 (1992). 

[23] G.I. Opat, S.J. Wark, and A. Cimmino, Appl. Phys. B 54, 396 (1992); T.M. 
Roach, H. Abele, M.G. Boshier, H.L. Grossman, K.P. Zetie and E.A. Hinds, 
Phys. Rev. Lett. 75, 629 (1995); A.I. Sidorov, R.J. McLean, G.I. Opat, 
W.J. Rowlands, D.C. Lau, J.E. Murphy, M. Walkiewicz and P. Hannaford, 
Quantum Semiclass. Opt. 8, 713 (1996). 

[24] S.J. Wark and G.I. Opat, J. Phys. B: At. Mol. Opt. Phys. 25, 4229 (1992); 
S.A. Schulz, H.L. Bethlem, J. van Veldhoven, J. Kuppcr, H. Conrad and 

G. Meijer, Phys. Rev. Lett. 93, 020406 (2004). 

[25] D. Cassettari, B. Hessmo, R. Folman, T. Maier and J. Schmiedmayer, Phys. 
Rev. Lett. 85, 5483 (2000). 

[26] E. Arimondo, H. Lew and T. Oka, Phys. Rev. Lett. 43, 753 (1979); P.E. 
Moscowitz, P.L. Gould, S.R. Atlas and D.E. Pritchard, Phys. Rev. Lett. 51, 
370 (1983); E. M. Rasel, M.K. Oberthalcr, H. Batclaan, J. Schmiedmayer 
and A. Zcilingcr, Phys. Rev. Lett. 75, 2633 (1995); D. M. Giltner, R.W. 
McGowan and S.A. Lee, Phys. Rev. Lett. 75, 2638 (1995); M. Kozuma, L. 
Deng, E.W. Hagley, J. Wen, R. Lutwak, K. Helmerson, S.L. Rolston and 
W.D. Phillips, Phys. Rev. Lett. 82, 871 (1999); S. Wu, Y.-J. Wang, Q. Diot 
and M. Prentiss, Phys. Rev. A 71, 043602 (2005). 

[27] Y.V. Baklanov, B.Y. Dubetsky and V.P. Chebotayev, Appl. Phys. A 9, 
171 (1976); Ch.J. Bordc, Ch. Salomon, S. Avrillier, A. Van Lerberghe, Ch. 
Breant, D. Bassi and G. Scoles, Phys. Rev. A 30, 1836 (1984). 

[28] M. Kasevich and S. Chu, Phys. Rev. Lett. 67, 181 (1991). 

[29] R.J. Cook and R.K. Hill, Opt. Comm. 43, 258 (1982); J.V. Hajnal and G.I. 
Opat, Opt. Comm. 71, 119 (1989); M. Christ, A. Scholz, M. Sciffer, R. 
Deutschmann and W. Ertmcr, Opt. Comm. 107, 211 (1994); A. Steanc, P. 
Szriftgiscr, P. Desbiolles and J. Dalibard, Phys. Rev. Lett. 74, 4972 (1995); 
R. Brouri, R. Asimov, M. Gorlicki, S. Fcron, J. Reinhardt, V. Lorent and 

H. Haberland, Opt. Comm. 124, 448 (1996); A. Landragin, G. Labeyrie, C. 
Henkel, R. Kaiser, N. Vansteenkiste, C. I. Westbrook and A. Aspect, Opt. 
Lett. 21, 1591 (1996); C. Henkel, K. Molmer, R. Kaiser, N. Vansteenkiste, 
C.I. Westbrook and A. Aspect, Phys. Rev. A 55, 1160 (1997). 



31 



[30] A.P. Kazantsev, Sov. Phys. JETP 40, 825 (1975); T. Slcator, T. Pfau, V. 
Balykin, O. Carnal and J. Mlynck, Phys. Rev. Lett. 68, 1996 (1992). 

[31] J. Oreg, F.T. Hioe and J. H. Ebcrly, Phys. Rev. A 29, 690 (1984); P. 
Marte, P. Zoller and J. L. Hall, Phys. Rev. A 44, R4418 (1991); P. Pillet, 
C. Valentin, R.-L. Yuan and J. Yu, Phys. Rev. A 48, 845 (1993); J. Lawall 
and M. Prentiss, Phys. Rev. Lett. 72, 993 (1994); P.D. Featonby, G.S. 
Summy, J.L. Martin, H. Wu, K.P. Zetie, C.J. Foot and K. Burnett, Phys. 
Rev. A 53, 373 (1996). 

[32] U. Gaubatz, P. Rudecki, M. Becker, S. Schicmann, M. Kulz and K. 
Bcrgmann, Chem. Phys. Lett. 149, 463 (1988); K. Bergmann, H. Theuer, 
and B.W. Shore, Rev. Mod. Phys. 70, 1003 (1998). 

[33] Y.B. Band, Phys. Rev. A 47, 4970 (1993); V.S. Malinovsky and P.R. 
Berman, Phys. Rev. A 68, 023610 (2003). 

[34] T. Pfau, Ch. Kurtsicfer, C.S. Adams, M. Sigel, and J. Mlynek, Phys. Rev. 
Lett. 71, 3427 (1993). 

[35] O. Houde, D. Kadio and L. Pruvost, Phys. Rev. Lett. 85, 5543 (2000); R. 
Dumke, T. Muther, M. Volk, W. Ertmer and G. Birkl, Phys. Rev. Lett. 89, 
220402 (2002). 

[36] B.W. Shore, K. Bergmann, J. Oreg and S. Rosenwaks, Phys. Rev. A 44, 
7442 (1991); K. Moler, D.S. Weiss, M. Kasevich and S. Chu, Phys. Rev. A 
45, 342 (1992). 

[37] Ch.J. Bordc, in Atom interferometry, edited by P. Berman (Academic, New 
York, 1997). 

[38] C. Cohen- Tannoudji, J. Dupont-Roc and G. Grynberg, Atom-Photon In- 
teractions: Basic Processes and Applications (Wiley, New York, 1992). 

[39] Ch.J. Bordc, in Fundamental Systems in Quantum Optics, edited by J. 
Dalibard, J.M. Raimond and J. Zinn-Justin (Elsevier, 1991); Ch.J. Borde, 
C. R. Acad. Sci. Paris, t. 2 (Serie IV), 509 (2001). 

[40] M. Kasevich and S. Chu, Appl. Phys. B 54, 321 (1992); B. Young, M. Kase- 
vich and S. Chu, in Atom interferometry, edited by P. Berman (Academic, 
New York, 1997). 

[41] N. Rosen and C. Zener, Phys. Rev. 40, 502 (1932). 

[42] Yu.N. Demkov and M. Kunike, Vest. Leningr. Univ. Fiz. Khim. 16, 39 
(1969). 

[43] S. Autlcr and C. Towncs, Phys. Rev. 100, 70 (1955); S. Gucrin and H.R. 
Jauslin, Adv. Chem. Phys. 125, 1 (2003). 

[44] T.-S. Ho and S.-I. Chu, J. Phys. B 17, 2101 (1984). 



32 



[45] U. Peskin and N. Moiseyev, J. Chem. Phys. 99, 4590 (1993). 

[46] V.S. Letokhov and V.G. Minogin, Zh. Eksp. Teor. Fiz. 74, 1318 (1978); 
Y. Castin and J. Dalibard, Europhys. Lett. 14, 761 (1991); M. Biichner, 
R. Delhuille, A. Miffre, C. Robilliard, J. Vigue and C. Champenois, Phys. 
Rev. A 68, 013607 (2003). 

[47] W. P. Schleich, Quantum Optics in Phase Space (Wiley, Berlin, 2001). 

[48] A. Iserles, H.Z. Munthe-Kaas, S.P. N0rsett and A. Zanna, Acta Numerica, 
215 (2000); R. Suarez and L. Saenz, J. Math. Phys. 42, 4582 (2001). 

[49] P.C. Moan and J.A. Oeo, J. Math. Phys. 42, 501 (2001). 

[50] X. Miao, Phys. Lett. A 271, 296 (2000). 

[51] M.V. Berry, Proc. R. Soc. London A 429, 61 (1990); K. Drese and M. 
Holthaus, Eur. Phys. J. D 5, 119 (1999). 

[52] Ch. Lubich, NIC Scries 10, 459 (2002). 

[53] G. Borrmann, Phys. Z. 42, 157 (1941); J. M. Cowley, Diffraction Physics 
(North- Holland, Amsterdam, 1990). 

[54] J. W. Knowlcs, Acta Cryst. 9, 61 (1956). 

[55] H. Rauch and D. Pctrascheck, in Neutron Diffraction (Springer, New York, 
1978). 

[56] Ch. Antoine and Ch.J. Bordc, Phys. Lett. A 306, 277 (2003). 

[57] M. Kasevich and S. Chu, Phys. Rev. Lett. 69, 1741 (1992); N. Davidson, 
H.J. Lee, M. Kasevich and S. Chu, Phys. Rev. Lett. 72, 3158 (1994). 

[58] CP. Agrawal, Applications of Nonlinear Fiber Optics (Academic Press, 
San Diego, 2001); Nonlinear Fiber Optics (Academic Press, San Diego, 
1995). 

[59] M. K. Oberthaler, R. Abfalterer, S. Bernct, C. Keller, J. Schmicdmayer 
and A. Zeilinger, Phys. Rev. A 60, 456 (1999). 



33 



